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Abstract 

A next-to-leading order QCD calculation for single-inclusive jet photoproduction in unpolarized 
^^^ • and longitudinally polarized lepton-hadron collisions is presented which consistently includes "di- 

pg , rect" and "resolved" photon contributions. The computation is performed within the "small-cone 

approximation" in a largely analytical form. Phenomenological aspects of jet production at future 
ep colliders such as the CERN-LHeC and the polarized BNL-eRHIC are discussed, placing partic- 
ular emphasis on the perturbative stability of the predictions and the possibility to constrain the 

'^i parton content of the photon. 

Oh' 
^ ; PACS numbers: 12.38.Bx,13.60.Hb,13. 
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I. INTRODUCTION 

Recent years have seen encouraging progress in the understanding of the inner structure 
of hadrons, triggered by precision measurements at high energy colhders, in particular the 
DESY-HERA. The knowledge gained on the dynamics which governs the interaction of 
color-charged particles is crucial as new experiments are probing even higher energies. The 
ultimate goal of the upcoming CERN-LHC, for instance, is the verification of the standard 
model by the discovery of a light Higgs boson or, else, the identification of new physics 
not anticipated within this framework. Subsequent measurements at a future linear collider 
will constrain the parameters of the scenario realized in nature with even higher accuracy. 
However, poorly understood QCD effects may render the interpretation of new signatures 
in terms of physics beyond the standard model difficult. The precise determination of the 
parton densities entering the description of any hadronic reaction is thus more timely than 
ever. 

Indeed, many improvements have been made recently on the parameterization of the 
proton distribution functions jj, [j, |3j. Much less is known about the hadronic structure of 
the photon. This ignorance severely limits the predictive power of photon-induced reactions 
since the resolved contributions, which are associated with its hadronic constituents, may 
be large at collider energies. As particularly promising means to cure this deficiency, jet- 
production processes in the photoproduction regime of lepton-hadron colliders, where the 
lepton beam acts as source of quasi-real photons, have been identified due to large production 
rates and small systematic uncertainties. The capability of HERA to constrain the partonic 
structure of the photon as well as the proton had been explored in great detail 4| and 
several NLO-QCD calculations for jet-photoproductionprocesses have become available, 

1 (sl, IgJ. For a polarization upgrade of 



which turned out to describe data reasonably we 

HERA, which has been discussed for some time |7|, a thorough sensitivity study has been 

jerformed in Ref. 8|. Polarized predictions at NLO-QCD accuracy have been presented in 

9j in the form of a fiexible Monte-Carlo program. 

Should single inclusive jets be used in the context of a "global analysis" including pho- 
toproduction data, however, fast codes are essential which are based on analytical methods. 
The aim of this work is thus to present a calculation of single-inclusive jet photoproduction 



in the framework of the "small-cone approximation" lOJ], which allows for an entirely ana 



lytical computation of partonic hard-scattering cross sections and has been demonstrated to 



approximate full jet cross sections extremely well in related pp-scattering reactions ll|, |l2| 



Predictions can then be made for photoproduction cross sections at future lepton-proton 



colliders such as the planned BNL-eRHIC |l3| and the CERN-LHeC |l4| which is cur- 



rently under scrutiny. In principle, the code developed also allows for the computation of 
jet-production observables at HERA. We refrain from presenting results for the HERA kine- 
matics here, however, since a vast number of NLO-QCD studies for this setting is available 

The plan of the article is as follows: In Sec. [Ill the technical framework used will be 
specified. Section IIIII contains numerical results for single-inclusive jet production at the 
LHeC and the eRHIC, respectively. Conclusions will be given in Sec. IIVI 

II. TECHNICAL FRAMEWORK 

We consider single-inclusive jet photoproduction in unpolarized and in longitudinally 
polarized lepton-proton collisions at NLO-QCD accuracy, i.e., the reaction ip -^ £'jetX. 
Provided the transverse momentum px of the jet is large, a polarized differential single- 
inclusive jet cross section can be written as a convolution 

= ^ / dXadXbAf^{Xa,fJ^f)Af^{Xb,IJ.f) 
a,b -^ 

X dAaab^}ctx{S,Xa,Xb,IJ,r,IJ'f), (l) 

where the subscripts "++" and "H — " refer to the helicities of the colliding leptons and 
protons, and S is the available c.m.s. energy squared. In Eq. ([T]), Xb denotes the momentum 
fraction of the proton which is taken by parton b and the Af^{xb, fif) are the longitudinally 
polarized parton distributions of the proton evaluated at a scale fif. The summation in 
Eq. ([T]) is performed over all partonic channels a + 6 — ;> jet + X contributing to the reaction 
ip — > £'jetX with the corresponding spin-dependent cross sections dA&ab^ietx, which are 
computed at NLO-QCD accuracy. 

The photoproduction cross section dAa which is observed in experiment consists of two 
pieces: First, the direct part dAadh, where a quasi- real photon emitted from the lepton 



beam scatters off parton b as an elementary particle such that a = 7 in Eq. ([T]). Second, the 
resolved contribution dAaj-cs, where the photon resolves into "hadronic" constituents, which 
in turn interact with the partons emerging from the proton. In dAaj-es, a denotes the parton 
stemming from the photon. Direct and resolved contributions can be cast into the form of 
Eq. dl]) by defining A/f as 

A/i(x„/.;) = r^ AP,,(2/) A/: U = ^,h) , (2) 
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+ 2m y^ [7^ --^2]) (3) 

VQmax rnjy^J J 

denoting the spin-dependent Weizsacker- Williams "equivalent-photon" spectrum for the 
emission of a circularly polarized coUinear photon with a virtuality smaller than Q^ax by a 
lepton of mass m^ J9J]. For the direct contribution, Xa has to be identified with the momentum 
fraction y of the lepton which is taken by the photon and thus 

A/; = 5{l - X,) . (4) 

In the resolved case, the A/J denote the parton distributions of the circularly polarized 
photon which are completely unmeasured so far. 

It is important to note that beyond the leading order neither (iAo"dir nor dAa-^cs are 
measurable cross sections per se, as their individual values depend on the factorization 
scheme chosen. Only if both pieces are evaluated using the same factorization prescription, 
their sum 

dAa = (iAcTdir + dAaj^cs (5) 

is a meaningful quantity. The unpolarized jet cross section, da = [da++ + da^_]/2, is 
obtained in complete analogy to the polarized one by replacing all spin-dependent parton 
distributions and partonic cross sections with their spin-averaged counterparts. The spin- 



averaged equivalent-photon spectrum can be found in Ref. |15[ |. 

In order to compute the hard-scattering cross sections d{A)aab^ietx, an algorithm has 
to be specified describing the formation of jets by the partons which undergo the hard 



scattering. A frequently adopted choice is to define a jet as the deposition of the total 
transverse energy of all final-state partons that fulfill 



iv - vT + 



b'? < R^ 



(6) 



where 77* and 0* denote the pseudo-rapidities and azimuthal angles of the particles and R 
the jet cone aperture. The jet variables are defined as 
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(7) 



Ej 



We resort to the so-called "small-cone approximation" [10|, which can be considered 
as an expansion of the jet cross section in terms of R of the form ^logi? + B + 0{R^). 
Neglecting 0{R^) pieces, the evaluation and phase-space integration of the partonic cross 



sections can be per 



has been shown 
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brmed analytically, as is our intention. The small-cone approximation 



12| to account extremely well for jet observables up to cone sizes of 



about R ^ 0.7 by explicit comparison to calculations that take R fully into account. The 
predictions of Ref. 12(] for single-inclusive jet production at the BNL-RHIC within the small- 
cone approximation are in excellent agreement with recent data from the STAR collaboration 

Since the resolved contribution (i(A)crres is technically equivalent to the jet production 
cross section in hadronic collisions, pp —>■ jetX, the corresponding partonic matrix elements 
squared can be taken from this previous calculation 121]. The direct contributions are 



adapted from the results for the d 



'A)d'^h^cX in single- inclusive hadron photoproduction 



171 . Il8l | with the techniques of Ref. 



12|. 



The hard-scattering cross sections are then implemented in a Monte-Carlo program which 
performs the convolutions with the parton distributions of the proton and with the equivalent 
photon spectrum numerically by means of an adaptive VEGAS integration. 

As consistency check, we have calculated jet cross sections for HERA kinematics and 
compared our results successfully to those of Refs. p, la, |9|. 



III. NUMERICAL RESULTS 

We now turn to a phenomenological study of jet production at future lepton-hadron 
colliders. Our aim is twofold: First, we will discuss the impact of NLO corrections on 
transverse momentum and rapidity distributions in various kinematic regimes and explore 
the stability of our predictions with respect to scale variations. Second, we will investigate 
the sensitivity of jet-production cross sections to the parton distribution functions of the 
photon and the proton and demonstrate how to increase the impact of specific contributions 
by adjusting parameters of the analysis. 

A. Single-inclusive jet production at the LHeC 

New opportunities for the measurement of single-inclusive jets in lepton-hadron colli- 
sions could be provided by a future LHeC ep collider at CERN, which is currently under 
scrutiny |14i] . In the following we will present predictions for the scenario where an electron 
beam circulates in the existing LHC tunnel with a nominal energy of Eg = 70 GeV, which 
in conjunction with the 7 TeV proton beam gives rise to ep collisions with a c.m.s. energy 
of v^ = 1.4 TeV. 

For the equivalent-photon approximation [l5| similar parameters as for the HI and ZEUS 
experiments at HERA are used, Q^^^ = 1 GeV^ and 0.2 < y < 0.85. All LO (NLO) 



calculations are performed with the CTEQ6L (CTEQ6M) parton distributions [19| and 
the according one-loop (two-loop) values for the strong coupling constant a^ in the MS- 
factorization scheme. For the parton distribution functions of the photon the LO (NLO) 



GRV set J20| is employed. 

The major motivation for studying jet cross sections beyond leading order is to reduce the 
theoretical uncertainties associated with a tree-level calculation. To access the improvement 
which can be gained by the inclusion of NLO corrections the single-inclusive jet cross section 
at the LHeC is shown in Fig. [1] as function of the jet transverse momentum px at LO and 
NLO. The jet cone size is set to i? = 0.7 and rapidities are integrated over —1 < r^ < 4 in the 
laboratory frame. The solid lines correspond to the setting /i^ = yUj = p^- The bands have 
been obtained by varying the factorization and renormalization scales simultaneously in the 
range pt/2 < fir = /i/ < 2pT- Throughout the pr interval considered, the scale dependence 
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FIG. 1: Single-inclusive jet cross section at the LHeC as function of px in NLO (multiplied by 
a factor of 10) and LO, integrated over —1 < rj < 4. The shaded bands correspond to a scale 
variation of the NLO and LO results, respectively, in the range pt/2 < ^r = /^/ ^ '^Pt- The lower 
panel shows the associated K factor. 

is very small at LO already, amounting to about 20 -r- 30 %. While at low values of pt 
the NLO corrections do not significantly improve the scale dependence of the cross section, 
towards higher pt the scale uncertainty of the NLO result is extremely small, going down to 
the level of 2.5 % aX pt = 200 GeV. This behavior indicates that the perturbative expansion 
is under excellent control, provided single-inclusive jets are produced at high transverse 
momentum. The impact of the NLO corrections on the cross section is indicated by the 
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FIG. 2: Relative contributions of different partonic subprocesses ab — > ietX to the NLO single- 
inclusive jet cross section at the LHeC, integrated over — 1 < r/ < 4. 



K factor, which we define as 



K(x) 



da^^^/dx 



(8) 



da^^/dx ■ 

For fir = fJ'f = Pt, K{pt) is larger than one everywhere, which refiects the increase of the 
cross section by the inclusion of the NLO contributions. 

Figure [2] illustrates the contributions of different partonic subprocesses to the rapidity- 
integrated NLO cross section. The resolved contributions are dominant over a large range 
of Pt with the direct contributions starting to take over only at pt ~ 200 GeV. Since the 
cross section is largest at low values of p^, this indicates that single-inclusive jet production 
at LHeC energies offers excellent opportunities for a more accurate determination of the 
parton content of the photon. 

To this end, the study of rapidity-differential cross sections is particularly suitable, since 
the momentum fractions of the hadronic constituents of the photon and the proton can be 
considered as functions of the rapidity of the observed jet in the laboratory frame, rj. As 
explained, e.g., in Ref. |8|, if counting positive rapidity in the forward direction of the proton. 
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FIG. 3: Relative contributions of different partonic subprocesses ab — > ietX to the NLO single- 
inclusive jet cross section at the LHeC, integrated over px > 20 GeV. 

large x-y — >■ 1 are probed at large negative values of 77. In this region, the direct contribution is 
expected to be largest and the f^ are dominated by the purely perturbative "pointlike" part 
which does not depend on the hadronic structure of the photon. The relative contributions 
of the various direct and resolved channels to the full cross section are depicted as function 
of rapidity in Fig. [31 

The LO and NLO rapidity-dependent jet cross sections integrated over pt > 20 GeV are 
shown in Fig. H] along with the associated scale uncertainties, obtained by varying renormal- 
ization and factorization scales simultaneously in the range Pt/2 < yUr = /^/ ^ "^Pt- Also 
shown are the relative contributions of the direct and the resolved cross sections at NLO. 
Strikingly, the scale uncertainty of the rapidity distribution is not improved by the inclusion 
of NLO corrections. This feature can be traced back to the large weight of contributions 
from relatively low pt, where the LO and the NLO scale dependences are of similar size, 
c.f. Fig. [H After imposing a transverse momentum cut of pt > 40 GeV, the scale depen- 
dence is slightly improved with the size of the cross section being reduced at the same time 
by approximately one order of magnitude. Even smaller scale uncertainties are obtained for 
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FIG. 4: Single- inclusive jet cross section at the LHeC as function of rj in NLO and LO, integrated 
over pt > 20 GeV. The shaded bands correspond to a scale variation in the range pt/2 < fJ-r = fJ-f 1^ 
2pT- The lower panel depicts the relative contributions of the direct and resolved sub-processes to 
the NLO cross section. 

Pt > 100 GeV, see Fig. [51 As the transverse momentum cut is increased, the impact of the 
direct photon contributions becomes larger and the resolved sub-processes start to be less 
important. 
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FIG. 5: Same as in Fig. HI but now integrated over px > 100 GeV. 



B. Polarized photoproduction of single-inclusive jets at eRHIC 



Photoproduction experiments complimentary to the measurements possible at the LHeC 
could be performed at the future lepton-proton collider eRHIC at BNL [13|. The polarized 
beams available at eRHIC offer unique opportunities for studying the spin structure of the 
circularly polarized photon, which is completely unknown so far. 

In the following we assume a hadronic c.m.s. energy of Vo = 100 GeV. For the equivalent- 
photon approximation [9|, [l5| we choose Q^ax — ^ GeV^ and 0.2 < y < 0.85. In the 
unpolarized case, we stick to the parton distributions of Sec. IIIIAI For the spin-dependent 
proton distribution functions we use the GRSV standard set [2V\ as default, since it contains 
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FIG. 6: Polarized and unpolarized single-inclusive jet cross sections at eRHIC as function of rj in 
NLO and LO, integrated over px > 4: GeV. The polarized cross section is multiplied by a factor 
of 30. 

both, a LO and an NLO parameterization. To illustrate the impact of different parton 
densities, we will also employ the new DSSV set [3'], which is only available at NLO, however. 
The parton distributions of the polarized photon are completely unmeasured so far. We 

22l | with minimal [A/'''(x, /iq) = 0] 



therefore consider the two extreme scenarios of Ref. 



and maxima 



[Af'^{x, Ho) = f^{x,fio)] saturation of the positivity constraint |A/'''(x,yUo)| < 
P^x^Ho) [23(1 . Here, /zq denotes the scale where the boundary conditions for the evolution 
are fixed. If not specified otherwise, the "maximal" set will be used. 

Figure [6] shows the rapidity dependent unpolarized and polarized single- inclusive jet cross 
sections at LO and NLO integrated over p^- > 4 GeV. The spin-averaged cross section 
receives sizeable positive NLO corrections resulting in 1.2 < K{ri) < L7. In the spin- 
dependent case, the NLO effects are rather small yielding a K factor close to one over a 
large range in rj. 

The scale uncertainty of the polarized cross section, illustrated by Fig. [TJ is sizably 
reduced at NLO. Scale dependences are generally smaller for jet cross sections than for 
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FIG. 7: Single- inclusive jet cross section at eRHIC as function of rj in NLO and LO, integrated over 
Pt > 4: GeV. The shaded bands correspond to a scale variation in the range prf^ < fJ-r = IJ-f < '^Vt- 

related hadron-production observables, such as the pion-production cross sections at eRHIC 
discussed in Ref . [l7| • This hierarchy is also observed in related pp-scattering reactions (see 



Refs. 



fl 



24j ) and is mainly due to the presence of scale dependent fragmentation functions 



in processes with identified hadrons. These non-perturbative objects describe the formation 
of hadrons from the final-state partons in a hard-scattering process at a specific scale, which 
may be different from the scale at which initial-state singularities are factorized into the 
distribution functions of the photon and proton, respectively. 

For extracting information on the hadronic structure of the circularly polarized photon 
the experimentally accessible spin asymmetry 



j^r. 



LL 



(9) 



da da^^ + da^ 

is most suitable. At large positive rapidities, A^^ is particularly sensitive to the parton 
content of the resolved photon as exemplified in Fig. [8], where the spin asymmetry is shown 
for the two extreme sets of polarized photon densities introduced above. To demonstrate 
that this sensitivity is not obscured by our currently incomplete knowledge of the polarized 
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FIG. 8: NLO-QCD spin asymmetry for single-inclusive jet photoproduction at eRHIC integrated 
over pt > 4: GeV for two different choices of proton distribution functions and the two extreme 
sets of polarized photon densities. 

proton, predictions are made for two different sets of proton distributions which mainly 
differ in the parameterization of the polarized gluon density. Towards negative values of t], 
All becomes less sensitive to the hadronic structure of the photon. With the qg subprocess 
being dominant in this region, the spin asymmetry could provide new information on the 
gluon density of the proton. 

IV. CONCLUSIONS 

In this work an NLO-QCD calculation of single-inclusive jet photoproduction in un- 
polarized and longitudinally polarized lepton-hadron collisions has been presented. The 
computation was performed in the context of the small-cone approximation. In this way, 
the evaluation of the partonic matrix elements and large parts of the phase-space integra- 
tion could be performed analytically, yielding a stable and much faster computer code than 
comparable programs that are based on a purely numerical approach. 
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We have performed a phenomenological analysis of jet photoproduction at future ep col- 
liders, in particular the LHeC and eRHIC We found that NLO corrections are sizeable at 
the LHeC with K factors of about 1.5. In the low-to-moderate pr regime the resolved photon 
contributions are by far dominant and scale uncertainties are large even beyond the leading 
order. If a large transverse momentum cut is imposed, these uncertainties can be brought 
down. At high p^, quark-initiated processes dominate over the gluonic channels. 

For eRHIC, we have focused on jet photoproduction by longitudinally polarized beams. 
Polarized cross sections exhibit smaller NLO corrections than their spin-averaged counter- 
parts with K factors close to one and only moderate scale dependences. The experimentally 
relevant double-spin asymmetry exhibits pronounced sensitivity to the hadronic structure of 
the (resolved) photon in the positive-rapidity regime. At negative values of 77, A^l could yield 
further information on the gluon polarization of the proton complementary to experiments 
in hadron-hadron collisions. 

In summary, future lepton-hadron colliders offer new possibilities for further constraining 
the hadronic structure of the (real) photon in the unpolarized case, and unique opportunities 
to access the completely unknown parton distributions of the polarized photon. Embedded 
in a global analysis of hadronic scattering reactions, single-inclusive jet production will 
also help to further pin down the parton densities of the proton, in particular its gluonic 
component. 
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